function fln  = fln(x,m,s)
fln=(1./(x.*s.*sqrt(2.*pi))).*exp(-(log(x)-m).^2./(2.*s.^2));
end

